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Abstract. The article deals with the mathematical model for media with hierarchical 
structure. Using the hamiltonian formalism, the dynamical system describing the state 
of hierarchically connected structural elements was derived. According to the analysis 
of the Poincare sections, we found the localized quasi-periodic and chaotic trajectories 
in the three-level hierarchical model. Moreover, studies of correlation functions shown 
that the power spectrum for three-level model possesses local maxima characterizing 
temporal scales with strong correlation. Using the Fourier analysis of the solution’s 
components, we have studied the distribution of energy injected in the system over 
hierarchical levels. Dynamical phenomena in the multi-level system were studied as 
well. 


1 Introduction 

One productive way to study dynamic processes in the lithosphere is an approach in which 
the lithosphere is considered as a hierarchical system of volumes (blocks) jIHS], where the 
block sizes are scaled in the interval from kilometers (tectonic plates) to millimeters (grains 
of rocks). The main features of dynamic behavior of such systems are high degree of nonlin¬ 
earity, complicated interaction between structural elements, the redistribution of energy be¬ 
tween hierarchical levels. These processes are manifested in the generation of high-frequency 
oscillations during the application of low-frequency disturbances in near - wellbore area jlj, 
in the change of spectra of aftershocks sequences [7j , in the existence of a wide range of non¬ 
linear waves after explosions in rock massifs [S] and so on. The exchange of energy between 
hierarchical levels is the significant factor in the loss of stability of block systems during 
seismic energy release mm- In order to study these dynamical phenomena, we propose a 
model for media consisting of the hierarchically connected elements (blocks). Most studies 
of the dynamics of hierarchical systems of oscillators have been focused on the synchroniza¬ 
tion processes [HH22], instead we are interested in exploring exchange of energy between 
hierarchical levels, in particular, energy transfer from the upper level to the lower one. 

The rest of this report is organized as follows. In Sec E] we introduce the models describing 
the hierarchically connected oscillators. Using the Poincare section technique, correlation 
and spectral analysis, in SecJ3] we investigate the simplified model for three-level system 
of oscillators. Sec. 0] is devoted to dynamical phenomena being observed in the multi-level 
hierarchical system when model parameters are varied. Concluding remarks are incorporated 
in final Sec El 
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2 Models for hierarchical media 


We model the hierarchical block medium by the embedded system of anharmonic oscillators 
(Fig. dK). The part of hierarchical system as a tree is shown in FigHJ). 

Assume that the model incorporates N levels and each oscillator from nth (1 < n < N ) 
level connects with s oscillators from the (n + l)th level, except the lowest level oscillators. 
The Hamiltonian of this system has the following form 
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where index (nk) denotes the oscillator’s position, namely, the nth level and the kth place; 
tk is the place of oscillator from the (n— l)th level connected with the (nk )th oscillator; m n k 
and p n k are the mass and, respectively, the momentum of the (nk )th oscillator; x n k is the 
oscillator’s coordinate; C n k is the stiffness of link, a is a constant. 
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Figure 1: The model of hierarchical medium (a) and corresponding tree (b). 
Inserting H into the equations 
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one can obtain the following system of differential equations 
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where sgn(x) = 1 if x > 0, otherwise sgn(x) = — 1. Considering the oscillators, which are 
identical on each level and move synchronously, system m can be written as follows 
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For convenience, let us introduce the new variables q n = x n — x n -\ having the meaning 
of displacements from the steady state, xq = const. Then x n = r A- Finally, system ((2]) 
reads 


^2<is = - <*>? | qif sgn (ft) + Iqi+if sgn(g m ). (4) 

a=i mi 

Using the designations F t = cof \q^ sgn(q , i ), </?* = sm^i/mj, we present system (J2|) in the 
form 

<?1 = — -^1 + 01 -^ 2 ; 

qn F'n—1 -Ui(1 T Vbi—l) T Fn-^n+li (5) 

QN = -PjV-1 _ Fn( 1 + <y?JV-l)- 

For simplicity, we assume that the quantities Ui and m, are the geometric sequences. So, 
uji = cuoU -1 , r > 1, and m; = mo/i* -1 , h < 1. From this it follows that (pi = sh = <p = const. 

3 Studies of the model with three hierarchical levels 

To begin with, let us consider the model with three hierarchical levels, i.e. putting N — 3. 
In this case system (EJ) has the following form 

qi = —Fi + p>F 2 , 

q 2 = Fi — F 2 (l + p>) + p>Fs, (6) 

^'3 = F 2 — F3 (1 + (p). 

Till now we did not achieve success in finding any exact solutions to system jBJ). Therefore, 
in order to get the information on the solutions structure, we are going to use the numerical 
and qualitative analysis methods. 

Due to the Hamiltonian nature of system ([HD, the trajectories of the dynamical system 
lay on constant energy hyper-surfaces. The allocation of these hyper-surfaces in the phase 
space can be effectively studied by means of the Poincare section technique [2j],27J. 

3.1 The sequence of frequencies = r z_1 , r = 1.1 

Let us fix the parameters ny = 1, r = 1.1, tp = 0.6 and choose the initial conditions for 
numerical integration in the form of q\ (0) = 0, q[( 0) = 0.2, r^-(0) = g'(0) = 0, j = 2,3. At 
first, consider the linear case when /3 = 1. Let the plane q\ = 0 be the Poincare section plane. 
Integrating of system ([3]) by means of the Dormand-Prince method [26j , we are interested in 
the points in which trajectories intersect the cross-plane in one direction only. Mapping this 
section onto the plane (q 2 ] qi), we obtain the diagrams presented in the figure Eh- The figure 
Efc> corresponds to the Poincare section of the trajectory beginning at the same starting point 
as in previous case but the parameter /3 = 1.05. 

In the same manner, putting f3 = 1.1, f3 = 1.15, /3 = 1.25, and /3 = 1.5 (the nonlinear 
cases), we obtain the typical Poincare diagrams depicted in the FigsEH31 
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Figure 2: Poincare sections at different (3 = 1 (a), (3 = 1.05 (b), (3 = 1.10 (c),/3 = 1.15 (d). 
Initial conditions q[ = 0.2 and the other coordinates are equal to zero. 


From the analysis of these figures it follows that the incorporation of nonlinearity in the 
model causes the appearance of the striped torus surfaces (FigsJ2b), tori with dense winding 
(FigsEb), trajectories similar to periodic ones (FigsJHb), which are not observed in the linear 
case. 

According to Figj3}l, the increasing of nonlinearity associated with the parameter (3 leads 
to the formation of Poincare sections uniformly filled with points. 

Comparing the initial conditions for construction of the Poincare diagrams of FigJ31 c and 
d, we see that small deviation of the first coordinate changes the observed regimes essentially. 
This is one of the main feature of the nonlinear models [28] • 

It worth noting that the sizes of regimes differ from each other weakly. But there are values 
of the parameter (3 when the region of the Poincare section (Fig®i,b,d) is filled with points 
nonuniformly. This tells us about existence of prevailing amplitudes in the quasi-periodic 
regime. 

Additional information on the properties of regimes we study can be obtain with the 
help of the correlation analysis. Using the solution of system (02) in the form of the discrete 
sequences q L (j ), i = 1, 2, 3, j = 1 • r, 2 • r..., M ■ r (r = 0.01 is the step of the temporal 
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Figure 3: Poincare sections at j3 — 1.20 (a), (3 = 1.25 (b), (3 = 1.5 (c), (d). Initial conditions 
for (a), (b), (c) q[ = 0.2 and other coordinates are equal to zero. Initial condition for (d): 
qi = 0.068, q[ = 0.2, q 2;3 = q 2 3 = 0 . In the corner of the figure (c) the quasi-periodic 
trajectory’s phase portrait corresponding to the Poincare section (c) is drown. 


variable discretization), we define the cross-correlation function [21] 
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If x = y, then we get the definition of the autocorrelation function R xx {k ). Together with 
the correlation functions we use their Fourier transformations FFT (/?.,,,,) = S xy (u>) which, in 
the case of autocorrelation function, coincide with the power spectrum of a signal [2fjJ. 

Let us begin from the correlation analysis of regimes derived at small (3. Taking the tra¬ 
jectories whose Poincare sections are depicted in FigsJ2fct,c, we obtain the graph of |,5j 1 (T)| 
(FigSK) and |,Sj .3 A) (Fig. 0b). According to Fig. 0b, the spectrum has the essential maxi¬ 
mum near u> = 0.1 and the other one at about 2u>. This tells us that the sequence q\ possesses 
the prevailing temporal scale when correlation is the most strong. Note that this dominant 
frequency displaces to the left at growing nonlinearity (to compare thick and thin curves). 
Since (3 is small, oscillators’ behavior does not differ from each other (Fig. 0b). 

For (3 = 1.5 when chaotic regime exists (Figj3ld), the spectrum of the function Rn has two 
multiple prevailing frequencies (Fig0b) which are accompanied by several peaks. At the same 
time, we observe more essential defragmentation of spectrum for the function i? i3 (Fig0]i). 
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Especially, it should be paid attention to the appearance of additional high-frequency peaks 
which can be regarded as supplementary temporal scales in media. 

3.2 The sequence of frequencies uji = r z_1 , r = 2.0 

Now we are interested in the studies of oscillations in system (131) when the partial frequencies 
are following a;* = r* _1 , r = 2.0. We keep the same initial conditions for numerical integration 
as in the previous subsection 10.1 1 

From the Poincare sections (Figl5h) it follows that the trajectories can form the pipe in the 
phase space at small values of parameter (3. When [3 increases, we observe the transformation 
of phase space in such a way that the regimes become more chaotic (Fig. [5b, c). Finally, the 
growth of (3 leads to essential chaotization of the Poincare section (Fig Si). 

The growth of (3 is thus accompanied by increasing of irregularity. In this situation we 
can not reduce system (J3J) to low dimensional model like one dimensional maps generated 
by Poincare sections. 

To understand phenomena in the system at increasing (3, we derive the Fourier transfor¬ 
mations for components q^, namely, FFT(qi ) at ft = 1.2 (FigEb-c) and FFT(qi ) at [3 = 1.5 
(FigJEb-f)- From spectra presented in FigEb-c it follows that all components q, are char¬ 
acterized by two essential spectral maxima. This causes the existence of torus (FigEb) in 
system’s phase space. 

For f3 = 1.5, the component qi has two wide spectral maxima (FigEb) and one notable 
maximum in the same place as in the spectrum from FigEb- 

The spectrum of qo (FigEb) involves three prevailing maxima having strongly irregular 
character. In the spectrum of q 3 we can distinguish three maxima as well. In contrast to the 
spectrum from FigEb, we see that the height of maxima grows as uj increases. This can be 
associated with the excitation of the high-frequency spectral modes. 

Therefore, system (El) can be regarded as a model describing the directed transport of 
energy in hierarchical media. 

4 Studies of the model with many hierarchical levels 

Consider now the model which consists of many hierarchical levels, for instance, N = 20. The 
fixed parameters are as follows uq = 1, r = 1.04, (3 = 1.15 and initial condition gi(0) = 0, 
q[( 0) = 0.3, qj{ 0) = g'(0) = 0, j = 2 ,..., N. Compare the qi and q^ components of solutions 
to system (El) derived at <p = 0.7 and p) = 0.9. Corresponding results are presented in the 
left and right panels of figures [7] and El respectively. 

For ip — 0.7 the amplitude of qi is much less than the amplitude of g iV • Snapshot from 
FigEb predicts that frequency anatomy of the qi and q\- signals is different. Moreover, the 
qi dynamics looks like a harmonic signal modulated by another signal with higher frequency. 
This is confirmed by the analysis of Fourier spectra for qi and q^ (FigEb)- The spectrum 
for qi depicted by thick line in FigEb contains one maximum corresponding to existence of 
mode like a cosine function and a number of much smaller maxima causing a ripple on the 
qi profile. 

As FigEb testifies, the spectrum for g iV depicted by the red thin curve contains a number 
of almost equivalent maxima located both in the low and high frequency zones. 
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Figure 4: The graphs of the functions |5ii(c(j)| (left) and |5i3(u;)| (right). The spectra 
correspond to the regime of Fig. Eb (top panels) and Fig. [31 (bottom panels). The curves 
depicted by thin lines correspond to linear system (EJ) with 3 = 1 . 
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Figure 5: Poincare sections at (3 = 1.2 (a), (3 = 1.3 (b), (3 = 1.4 (c), /3 = 1.5 (d). The values 
of parameters r = 2, ip = 0.6. Initial conditions q[ = 0.2 and the other coordinates are equal 
to zero. 
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Figure 6: (color online) Fourier spectra of the q L , i = 1,2,3-component of the trajectories 
drawn in the Fig. [5] at (3 = 1.2 and 3 = 1.5, respectively. 

Similar analysis can be carried out for the solutions to system (J3J) at p = 0.9. Snapshots for 
the components qi and q^ plotted in Fig0) show that the amplitude of <ji is larger than q^. 
According to FigEb, the main part of q\ Fourier spectrum is localized in low frequency zone. 
But the spectrum of q n is almost uniformly distributed in the domain under consideration. 

From the results presented in this section we can conclude that, at first, a critical value 
of tp can exist corresponding to the formation of comparable oscillations on the first and 
the last hierarchical levels of media. Secondly, the spectrum of the lowest level we are most 
interested in is distributed in a wide frequency domain and depending on the p dominant 
frequencies can be distinguished. 

5 Conclusions 

In summary, the strongly nonlinear system of coupled oscillators describing media with 
hierarchical structure is introduced. The problem we have stated is concerned with the 
peculiarities of energy transfer along the hierarchical structure of media. To treat this 
problem, we have studied the phase portraits, correlation functions and Fourier spectra 
which characterize oscillators placed on different hierarchical levels of media. We thus have 
seen that hierarchical systems manifest quasiperiodic and chaotic regimes development of 
which depends on the auxiliary parameter ip. It is shown that among values of the parameter 
ip corresponding to processes of transferring energy from the top level to the lowest one there 
is a threshold value. Therefore, within the framework of presented model, one can confirm 
that the hierarchical structure accompanied by nonlinearity plays an important role in the 
transformation of energy flows in media. Rearrangement of the hierarchical levels caused 
by intensive loading is the natural mechanism of accumulation and radiation of energy in 


9 























































Figure 7: (color online) Snapshots of the components qi(t) (thick curve) and (pv(t) (thin red 
curve), t e (330; 500) at <p = 0.7 (a) and p = 0.9 (b). 
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Figure 8: (color online) Fourier spectra for the components q\(t ) (thick curve) and quit) 
(thin red curve) at <p = 0.7 (a) and p = 0.9 (b). 
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structured media under strongly non-equilibrium conditions. 
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